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^ ' Abstract^: This paper reproves a few results concerning paperfolding 

sequences using properties of Catalan numbers modulo 2. The guiding prin- 
C^^ , ciple can be described by: Paperfolding = Catalan modulo 2 + signs given 

by 2— automatic sequences. 
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1 Main results 

In this paper, a continued fraction is a formal expression 

^0+ ^ , ^2 = ^0 + ai/(bi + 02 ih2 + . 
"1 + 62+. •• 



^ ■ with 6o; ^ii • • • 5 oi, 02, . . . G C[x] two sequences of complex rational functions. 

O . We denote such a continued fraction by 

^D ' and call it convergent if the coefficients of the (formal) Laurent-series ex- 

c~| ■ pansions of its partial convergents 

a 1^0 _ Oo n _ OoOi+oi ^k _, Oil a2| , , a/ci 

Qo 1 Qi hi Qk \bi I62 \ok 

> 

K> ', are ultimately constant. The limit "^k^n^k^^ £ C![[x]][x~-'^] is then the 

5-H I obvious formal Laurent power series whose first coefficients agree with the 

Laurent series expansion of Pn/Qn for n huge enough. 

Let Wi be the word {—xi) xi of length 2 in the alphabet {xi, —xi}. For 
k >2, we define recursively a word Wk of length 2(2^^' — 1) over the alphabet 
{ibxi, ±X2, . . . , ±2;^} by considering the word 



Wk = Wk-i Xk i-Xk) Wk-i 



^Math. Class; 11B65, 11B85, 11C20, 15A23. Keywords: Catalan numbers, paperfold- 
ing, automatic sequence, continued fraction, Jacobi fraction, binomial coefficient. 



constructed by concatenation where Wk-i denotes the reverse of Wk^i ob- 
tained by reading W^-i backwards. The first instances of the words Wk 
are 

Wi = {-xi) xi 

W2 = {-Xl) Xi X2 {-X2) Xi {-Xl) 

W^3 = {-Xl) Xl X2 {-X2) Xl {-Xl) X-i {-X-i) {-Xl) Xl {-X2) X2 Xl {-Xl) 

and these words are initial subwords of a unique, well-defined infinite word 
Woo = wi W2 W3, . . . with letters Wi G {±xi, ±X2, • • •}• 

The following result has been known for some time, cf [S] or J13j . Similar 
or equivalent equivalent statements are for instance given in [0], Theorem 
6.5.6 in P]. It is closely related to paperfolding, see [2], JJ, [3]. Related 
results are contained in [3j, ^ and [T!^ . 

Theorem 1.1 Consider a sequence xi,X2, ■ ■ ■ with values in C If the cor- 
responding continued fraction 

ii _L !fli + !f^ + I ^fcl I 
|i |i |i ■■■ |i 

associated with the sequence Woo = wi W2 ■ ■ ■ defined by the letters of the 
infinite word Woo converges, then its limit is given by 

00 j 
H- Xl -I- x1x2 + xfx^xs + xf X2X3X4 -I- . . . = 1 -I- 2J Jl 2;f ' . 

Remarks, (i) In fact, the proof will follow from the identity 



2-'- 



^2(2^-1) 11 11 11 11 Pit=l 

(cf. Lemma 6.5.5 in 2 ) showing that the 2(2 — 1)— th partial convergent 
of this continued fraction is a polynomial in xi, . . . , x^. 

(ii) It is of course possible to replace the field C by any other reasonable 
field or ring. 

Examples. (1) Setting Xj = x for all i and multiplying by x we get 

nf* rf* /y* /y /y» ry* 

*Xj Ju ti-/ t-ly <Xy tXj 

00 

2 4-8 X ^ 2^ 

= x-|-x + x + X -|-... = /^x 

fc=0 

The sequence of signs wi,W2, ■ ■ ■ € {±1} given by 

00 



X = —+Wi—+W2l-+W^— + 



can be recursively defined by 

W8i+3 = -W8i+A = (-1)' , 
■"^82+7 = -Wsi+S = W4i+3 

and is 2— automatic (cf. 2 ). 

(2) Setting Xi = x^^^ and multiplying by x we get 

I 2l 2l 4l 4l 2l 

tXj tXi tXi tXj tiy 



= X + X^ + X^ + X^^ + . . . = ^ X^ . 

fc=0 

(3) Setting Xj = x-^+(*~-^)*- and multiplying by x we get 

I I I Si Si I 

ity Jy Jy ti/ t*/ 

|T~|T'^|T'^li ir^Ji^'" 



X + x^ + x*^ + x^^ + ■ . . = y^ 



fc! 



X 
k=l 

Consider the sequence fj, = (fiQ, fj,i, ■ ■ ■) = (1, 1,0, 1,0,0,0, 1, . . .) defined 
by ^Y^'^Q^ix'^ = YlT=o^'^ ■ ^^ associate to /i = (/io, /Ui, • • •) its Hankel 
matrix H whose (i, j)— th coefficient /ijj = Hi+j-, < i,j depends only on 
the sum i+j of its indices and is given by (i+j)— th element of the sequence 

Consider the (2— automatic) sequence s = (sq, si, S2, • • •) = 1 1 (~1) 1 (~1) (~1) 
defined recursively by sq = 1 and 

S2i = (-1)* Si , 
S2i+1 = Si . 

It is easy to check that s„ = (— l)^o(") where Bq^u) counts the number of 
bounded blocks of consecutive O's of the binary integer n (or, equivalently, 
Bo(n) equals the number of blocks 10 appearing in the binary expansion 
of the n). E.g, 720 = 2^ + 2^ + 2^ + 2^ corresponds to the binary integer 
1011010000. Hence -Bo(720) = 3. This description shows that we have 
•S2"+i+n = ~Sn ^ov n < 2". Since we have Bq{2'^~^^ + n) = 1 + BQ{n) for 
n < 2" and So(2"+i + n) = Bo{n) for 2" < n < 2"+\ the sequence sq, si, . . . 
can also be constructed by iterating the application 

where where Wa, W^ are the first and the second half of the finite word 

So . . . S2"-l- 



Denote by D3 the infinite diagonal matrix with diagonal entries sq, si, S2; • 
and by Da the infinite diagonal matrix with diagonal entries given by the 
alternating, 2— periodic sequence 1, —1, 1, —1, 1, —1. 

We introduce the infinite unipotent lower triangular matrix L with co- 
efficients ^ij G {0, 1} for < z, J given by 



«j 



2i 
I- J 



2i 
- 3 - 



2i 
i- 3 



+ 



2i 
■ 3 - 



2i + l 



(mod 2) 



L(8) 



and denote by L{k) the k x k submatrix with coefficients lij, < i,j < k 
of L. The matrix L(8) for instance is given by 

/I \ 

1 1 

Oil 

1111 

11 

1111 

110 11 

V 1 1 1 1 1 1 1 1 y 

Theorem 1.2 (i) We have 

H = Ds L Da L' Ds . 

(a) (BAOBAB— construction for L) We have for n > 

\ 

D A 

Lf2"+2^ 



.(2"+^) 



1 ^ 

B 


A 







B 

A 


A 
B 



A J 



where A and B are the obvious subniatrices of size 2" x 2" in L(2""'"^). 

Remarks, (i) Assertion (ii) of Theorem 1 1 . 21 vields an iterative construc- 
tion for L by starting with L{2) = I 

(ii) Conjugating H with an appropriate diagonal ±1— matrix, we can get 
the the LU-decomposition of any Hankel matrix associated with a sequence 
having ordinary generating function x~^ [Ylk^o ^^x^ ) ' ^0 = !> ^i) £2, • • • £ 
{±1}. More precisely, the appropriate conjugating diagonal matrix has co- 
efficients dn^n = [—ly^3=o^j'^i+i for < n = '^j=Qi'j2^ a binary integer. 
The case eo = — 1 can be handled similarly by replacing the diagonal matrix 
Da with its opposite —Da- 

The matrix P = (Dg L Dg)^ = Dg L^^ Dg with coefficients pi^j, < i,j 
encodes the formal orthogonal polynomials Qo,Qi, ■ ■ ■ ,Qn = Ylk=oPn,kx'' 



with moments /^o, /^i, • • • given by ^ 



00 i+l 

f^iX 



E 



k=0- 



In order to describe the matrix P we introduce a unipotent lower trian- 
gular matrix M with coefficients mjj- € {0, 1}, < i, j by setting 



rrii^j = f 2 • ) (™°"^ 2) • 



As above (with the matrix L) we denote by M{k) the k xk submatrix with 
coefficients w-jj-, < i, j < A; of M. 

Theorem 1.3 (i) We have 

LDaM = Da 

(where Da denotes the diagonal matrix with 2— periodic diagonal coefficients 

1,-1,1,-1,...;. 

(a) (BAOBAB— construction for M) We have have for n > 



M(2"+^) 



(A' \ 

B' A' 

A' B' A' 

\ B' B' A' I 



r^+l^ 



where A' and B' are the obvious submatrices of size 2" x 2^ in M{2 

The matrix P = Ds Da M Da Ds = (D^ L Ds)'^ encoding the coeffi- 
cients of the orthogonal polynomials with moments /xo,/ii, . . . has thus all 
its coefficients in {0, ±1}, see also P and [S]- 

The (formal) orthogonal polynomials Qo,Qi, ■ ■ ■ satisfy a classical three- 
term recursion relation with coefficients given by the continued Jacobi frac- 
tion expansion of ^ ^Uj x* = ^ ( Y^'^=o ^^ ) • This expansion is described by 
the following result. 

Theorem 1.4 We have 

oo 

E2'' _ x\ x^l , x'^\ 

•^ ~ \l-dix "•" \l-d2x "•" \l-d3x "!"••• 
A;=0 

-^1 I -^ I I ^ \_ I ^£ \_ I -^ I I 

~ |l-x "•" |l+2a; "•" |1 "•" |1 "•" |l-2x "•" • • • 

where dn = ^"~^"~"^ with s_i = 0, sn = 1 and S2i+i = (— l)*S2i = Si is the 
recursively defined sequence introduced previously (extended by S-i = 0). 

The formal orthogonal polynomials with moments ^Q,fj,i,... are thus 
recursively defined by(5o = li Qi = x + 1 and 

Qn+i = (x- dn)Qn + Qn~i for n > 2 . 



The rest of the paper is organised as fohows. 

Section 2 contains a proof of Theorem 11.11 

Section 3 recaUs a few definitions and facts concerning Hankel matrices, 
formal orthogonal polynomials etc. 

Section 4 is a little digression about Catalan numbers. All objects ap- 
pearing in this paper are essentially obtained by considering the correspond- 
ing objects associated with Catalan numbers and by reducing them modulo 
2 using sign conventions prescribed by recursively defined sequences (which 
are 2— automatic). 

Section 5 recalls two classical and useful results on the reduction modulo 
2 (or modulo a prime number) of binomial coefficients. 

Section 6 is devoted to the proof of Theorem 11.21 

Section 7 contains a proof of Theorem 1 1 . 31 and a related amusing result. 

Section 8 contains formulae for the matrix products M L and L M. 

Section 9 contains a proof of Theorem 11.41 

Section 10 displays the LC/— decomposition of the Hankel matrix H as- 
sociated with the sequences 1,0,1,0,0,0,1,0,... obtained by shifting the 
sequence 1, 1, 0, 1, 0, 0, 0, 1, . . . associated with powers of 2 by one (or equiv- 
alently, by removing the first term). 

Section 11 contains a uniqueness result. 

2 Proof of Theorem 11.11 

More or less equivalent reformulations of Theorem 11.11 can be found in 
several places. A few are jH], |llj and [2]. We give her a computational 
proof: Symmetry-arguments of continued fractions are replaced by polyno- 
mial identities. 

Given a continued fraction 

, ail ^2! 

it is classical (and easy, cf. for instance pages 4,5 of ^01) to show that the 

product 

1 bo\ f ai\ f a2\ ( at 
1 M 1 bi J [1 b2 J"'[l bk 



equals 



where 



Pk-i Pk 

Qk-l Qk 



Pk , Ol , , a-kl 

Qk \bi \bk 



denotes the fc— th partial convergent. 



Given a finite word U = ui U2 ... ui we consider the product 

Ui \ f U2 \ f Ul 



^(^)= VI 1 M 1 1 / VI 1 
Defining 

Xk = i + EUUU 

Xk = i+E-=im=i-f"-nti^^ 

(example: X3 = 1 + xi + x\x2 + xf X2X3 and X3 = 1 + xi + xf X2 — x^x^xs) 
one checks easily the identities 

Xn + Xn+l{X^_]^ — XnXn) = Xn+l 

and 





Xn — Xn+l{X^_i — XnXn) — Xn+l ■ 


Lemma 2.1 


We have for all n > 1 




V -^n-1 ^n J 


and 









Proof. The computations 

— xi \ / xi \ / — xi — xi 



M(iy,) = M((-x,)xi; > ^ 1 Ml 1 ; VI i+x, 



(l-Xl)-l2 1-(1 + Xl] 

and 



1^ 1 + Xl 



Xl \ / -Xl \ _ / (l+xi)-l2 l_(l-xi) 

1 1 J V 1 1 J ~ V 1^ l-2;i 



M{Wi) = M(xi (-Xl)) = 

prove Lemma I2. II for n = 1. 
For n > 1 we have 

M{Wn+l) = M{Wn Xn+l (-X„+l) W^) 

^n — ^n-l 1 ~ -^n \ f ^n+1 Xn+l \ ( Xn — -^„_i 1 — Xn 

xi^i x„ ;v 1 i-x„+i;v K-i Xn 

{Xn — Xn+l{Xn^i — XnXn)) — Xn 1 — {Xn + Xn+l{X^_-^ — XnXn)) 

Xn Xn + Xn+l{X^_'^ — XnXn) 

Xn+l — Xn 1 — Xn+l 
Xn Xn+l 



and 



M{Wn+l) = M{Wn (-X^+i) X„+i Wn) 

Xn — ^n-1 1 ~ Xn \ ( —Xn+l —Xn+1 \ ( Xn — ^„_i 1 — X„ 



Xtl Xn J\ 1 1 + X„+i ; V X2_, Xn 

{Xn + iCn+l(^n-l ~ -'^n-'^^n)) " -''^n ^ ~ (^n ~ ^n+l(-^n-l ~ XnXn)) 

Xn Xn — Xn+l{X^_-^ — XnXn) 

X„M1 — X„ 1 — Xr 



X„ X 



n+1 



which proves Lemma l2. II bv induction on n. D 

Proof of Theorem 11.11 Lemma I^TTl shows that 

M{lWn) = 



1 ^ 


\( Xn-Xl_, 


1-X, 


1 ly 


' V xl_. 


Xn 




Xn \ 




Xn 


1 J 





The 2(2"" — 1)— th partial convergent of the continued fraction 

1| Wl\ W2\ 

|T "^ IT "^ IT "^ ■ ■ ■ 

defined by the word W^o equals thus 

n j 



j=i k=i 



23 -k 

^k 



and the sequence of these partial convergents has limit 

oo j 

J=l fc=l 

which is the statement of Theorem 11.11 D 

3 Formal orthogonal polynomials 

The content of this section is classical and can for instance be found in |3j 

or [in]. 

We consider a formal power series 

oo 

fc=0 
as a linear form on the complex vector space C[x]* of polynomials by setting 

n n 

i=0 i=0 



We suppose that the n— th Hankel matrix H{n) having coefficients hij = 
f^i+j, < i,j < n associated with ^ = (^O)---) is invertible for all n. 
Replacing such a sequence ji = (//q, A*i, • • •) by -^/U = (1, ^, . . .) we may 
also suppose that ^q = 1. We get a non-degenerate scalar product on C[x] 
by setting 

{P,Q)=lg{PQ). 

Applying the familiar Gram-Schmitt orthogonalisation we obtain a sequence 
Qo = 1, Qi = X - ao, 

Qn+l = {x - an)Qn " KQn-l, n > 1 

of mutually orthogonal monic polynomials, called formal orthogonal polyno- 
mials with moments fio, ^i, . . .. 

We encode the above recursion relation for the polynomials Qo,Qi, ■ ■ ■ 
by the Stieltjes matrix 



S 



( ao 1 
hi ai 1 

V 



\ 



/ 



The non-degeneracy condition on hq, fii, . . . shows that the (infinite) 
Hankel matrix H with coefficients hij = iJ,i+j, < i,j has an L{7— decomposition 
H = LU where L is lower triangular unipotent and U = DL^ is upper tri- 
angular with D diagonal and invertible. 

The diagonal entries of D are given by 



and we have thus 



fc=i 



n— 1 



det{H{n)) = ll{bkr 



-k 



k=l 



where H{n) is the n— th Hankel matrix with coefficients hij 
i,j < n. The i— th row vector 

(mi,o,"^j,l,...) 



fJ'i+j, < 



of M = L ^ encodes the coefficients of the i—th orthogonal polynomial 

i 

Qi = '^nii^kx'' . 



fc=o 



The Stieltjes matrix S satisfies 



LS = L^ 



hfi h,i h,2 



where L_ is obtained by removing the first row of the unipotent lower trian- 
gular matrix L with coefficients /jj-, < i, j involved in the LL'^— decomposition 
oiH. 

The generating function g = X^fcLo ^^kx^ can be expressed as a continued 
fraction of Jacobi type 

g{x) = 1/(1 - aox - bix^/{l - aix - b2x'^/{l - a2X - 63x^/(1 - . . . )))) 

where oq, ai, • • • and 61, 62, • • • are as above and are encoded in the Stieltjes 
matrix. 



4 The Catalan numbers 

The Catalan numbers C*™ = ( ^) ^q::i are the coefficients of the algebraic 
generating function 



C{x) = En=oCnX'' = ^'^^ = 1 + X + 2x^ + 5x^ + Ux^ + ■ , 

1| x\ x\ x\ 

|T~]T~]T~]T~'" 

11 x^l x^l x^l 



\l-x |l-2x |l-2x |l-2x 

and satisfies c = 1 + x c^ . 

We have g{x) = c{s) (mod 2) where g{x) = YlTLo^"^^'^^ '^^- ^^^ ^^^^ 
lines of ^J! . A different proof can be given by remarking that Cn counts the 
number of planar binary rooted trees with n + 1 leaves. Call two such trees 
equivalent if they are equivalent as abstract (non-planar) rooted trees. The 
cardinality of each equivalence class is then a power of two. An equivalence 
classe containing only one element is given by a rooted 2— regular trees 
having 2" leaves which are all at the same distance n from the root. 

Many interesting mathematical objects associated with g{x) are in fact 
obtained by reducing the corresponding objects of c(x) modulo 2. Mirac- 
ulously, all this works over Z by choosing suitable signs given by a few 
recursively-defined sequences (which are in fact always 2— automatic, see j2] 
for a definition). 



10 



Consider the unipotent lower triangular matrices 



L 



and 



L 



/ 1 








\ 


1 


1 








2 


3 


1 






5 


9 


5 


1 




14 


28 


20 


7 


1 


^42 


90 


75 


35 


9 l) 


1 










2 


1 








5 


4 


1 






14 


14 


6 


1 




42 


48 


27 


8 


1 



V 132 165 10 44 10 1 / 

^^. - u'-,) - G4-i)' respectively k, = i^l^j) - G^^J 
< i,i obtained by taking all even (respectively odd) lines (starting with 
line number zero) of the so-called Catalan-triangle 



with coefficients h 



for 



1 

1 1 

2 1 

2 3 1 

5 4 1 

5 9 5 1 

14 14 6 1 

14 28 20 7 1 



The matrix L is associated with the LC/— decomposition of the Hankel matrix 
H = L L^ having coefficients hij = d+j = i^^^p)/{i + j + 1), < i,j 
the Catalan numbers. The matrix L yields the LC/— decomposition of the 
Hankel matrix H = L U having coefficients hij = Cj+j+i = ( It'^ti )/{i + 
j + 2), < i,j associated with the shifted Catalan sequence 1, 2, 5, 14, 42, .... 
The inverse matrices L~^ and L~^ are given by Da M Da and Da M Da 
where Da is the diagonal matrix with 2-periodic alternating diagonal coef- 



ficients 1, —1, 1, —1, and where M and M have coefficients m 



*j 



m- 



\ 2j+l ) 



m 



and 



M - \2j+il for 0<i,j. 
Let us mention that the products P 



LM and P = L M have non-zero 



coefficients pij 



(2»)! j! 



i\ (2j)\ (i-jy. 



and pij = 4^* ^^ (*) for < j < i. They can 



11 



also be given by 



P = L M = exp( 



2 
6 
10 

V 



) 



and 



P = L M = exp( 



4 



V 








12 



■•• / 



Computations suggest that the products M L and M L have also nice 
logarithms: 



ML 



exp( 



2 

6 

2 10 

6 

2 10 

6 

V : : 



\ 





14 



14 




18 





22 



) 



and 



M L = exp( 



/ 










4 














8 









4 





12 










8 





16 





4 





12 





20 






V : 



5 Binomial coefficients modulo 2 

Chercher ref. (Concrete Maths? Graham, Knuth, Patashnik) 
Using the Frobenius automorphism we have 

(l+x)^^"^^^ =Y[{l+x'^^y^ (mod 2) 



12 



which shows the formula 



).n(:) (-^2, 



where n = Yli>o ^i 2*' ^« ^ {*-*' ^^ ^"^^ ^ ~ Z]i>o ^« ^'' '^i ^ i^' -'^i ^^^ ^^^ 
natural binary integers. Kummer showed that the 2— valuation V2 of (^) 
(defined by 2^2|(^) and 2^^+^ /f(fc)) equals the number of carry-overs when 
adding the two natural numbers k and n — k, written in base 2 (cf. 16~). 
As an application, one (re)proves easily that Cn = {^) ;^ = 1 (mod 2) 

ifandonlyifn+1 G {2^}fceN is apower of two. Indeed, (^^);^ = nf'^n+iy. = 
{^^ ) (mod 2) and addition of the binary integers n and {n + 1) needs 
always a carry except if n = Yli=o '^ ^^^ n+ I = 2^. 

6 Proof of Theorem 11.21 

The aim of this section is to prove Theorem ll. 21 which describes the LC/— de- 
composition of the Hankel matrix H with coefficients hij = Hi+j, < i,j 
where 

oo oo 

^/ijX* = ^2;^'^^ = l + x + x^ + x'^ + x^^ + .... 

k=0 j=0 

We give in fact 3/2 proofs: We first prove assertion (i) by using the 
definition kj = (^i^^) (mod 2) of L. 

We prove then assertion (ii) showing that the matrix L satisfies the 
SAO-B^i?— construction. 

Finally, we reprove assertion (i) using the recursive i?AOi?yli?— structure 
of L. 

Proof of assertion (i) in Theorem II. 2L We denote by (^) G {0, 1} 
the reduction (mod 2) with value in {0, 1} of the binomial coefficient (Tj. 
The main ingredient of the proof is the fact that 



'2n + 1\ _ / n 

k ) = [lk/2\^ 

where [fc/2j = k/2 if k is even and [k/2\ = {k - l)/2 if k is odd. 

We compute first the coefficient r2i^2j of the product R = DgLDaL^Dg 
(recall that L has coefficients kj = (/.t)): 



r2^,2, = ^2..2,E.(-l)^(|gKSSl ^^^^_ 

= S2iS2j ( Z^fc=0 \2i-2k) (2J-2J ~ Z^k=l \2i-2k+l) {2j-2k+l) ) 

= S2iS2j {J2k=0 [i-k) [j-k) ~ ^k=l \i-k) \j-k)) 
_ (2i\ (2j\ _ / 1 if i = j = 

- S2^S2j{^){^) -| Q otherwise. 

13 



Similarly, we get 



fc(.-i j bi+l-fcj bj+l-fcj 

= S2J+lS2i+l ( 2L/fc=0 (2i+l-2j (2J+I-2J ~ /^fc=0 (2i-2fc) (2j-2fc. 

For r2i,2i+i = r2j+i,2i we have 

/v^ / 4i+l \ / 4 1+3 7 v^ 7 4i+l \ / 4i+3 \\ 
r2i,2i+l - g22g2j+l ^2L,fc \2i-2k)\2j+\ ~2k) ~ Z^fc l2»--2fc-lj l2j-2fcjj 

= S2iS2j+l ( Efc (i!fe) (^^fc ) - Efc G-fc-l) (^/^fc 



S2iS2i+l (Efc (G-fc) + G- fc-l) j ( A) -^Efc (i-fc-l)(j^fc 
/^^ ('2i+l\ /2j+l\ r, v^ /' 2i \ /2j+l 

Using parity arguments and the recm'sion s^i = (— l)*Si, S2j+i = Sj we get 
for i even 



r2,2,+l = ^2.^2,+! E.(-l)'^(^^)(^) + 



,2V ( 2*+^ ) - ( 2» ]] ( 2j+i 

T^^Z^fc Ui-2fc-ly' \i-2k-2j I \j-2k-l 



Similarly, for i odd we have 



r2,2,+l = -«2..2,+l E.(-l)'=(^!^)(^) + 



r, V- / /2i+l\ _ / 2i \ \ /2j+r 
^ l^k \ \i~2k) \i-2k-l) j \j-2k. 



-.-. E.(-i)^(^!^)(^) + o 



j~k) ^ ^ ) — 'hi ■ 

Finally, the above computations together with induction on i + j show that 

rjj = except if i + j = 2^^ — 1 for some fc G N. In this case we get rjj = 1 

which proves the result. □ 

Proof of assertion (ii) in Theorem II. 2L Assume that L(2"+^) is 

given by I . \. (This holds clearly for n = 0.) 

We have for < i, j < 2"+^ 

since 2i + 1 < 2""*"^. This shows already that L(2"^^) is of the form 

^ J 1 with A = L(2"+i^ 
B A 
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Consider now a binomial coefficient of tlie form 



2(2"+i + i) + 1 
(2"+i +i) -j 



2"+2 + 2i + l 
2«+i + i- j 



with < i,j < 2"+ which is odd. This implies i + j + 1 > 2"+ since 
otherwise (2"+i + i)-j < (2"+i+i)+j + l < 2"+^ < 2"+^ + 2i + 1 and there 
must therefore be a carry when adding the binary integers ((2"+^ + ~ j) 
and ((2"+i + i)+j + l). 

Two such binary integers ((2""+^ + ^) — i) and ((2"+^ + i) + j + 1) add 
thus without carry if and only if the binary integers (2"+^ + 2""''^ + i 



and (2"+^ + i + j + 1 



in+2\ 



add without carry. This shows the equality 



2(2 



ra+l 



+ i) + l 



(2"+i + 



J 



2(2"+i+i) + l 
;2«+i + i) - (2"+2 - 1 



J, 



for < i,j < 2"^"'^. Geometrically, this amounts to the fact that the block 
B is the vertical mirror of the block A and this symmetry is preserved by 
the Byl05^5-construction. D 

BAOBAB— proof of assertion (i). We have 



1 1 
1 1 

V 1 1 1 1 / V 



■1 / 



\ /I 1 1\ 

1 1 1 
1 1 



V 



/ 1 1 1\ 
10-10 
0-100 

V 1 0/ 



and conjugation by 
/I 



v 



yields 



1 / 



/ 1 1 1\ 
10 10 
10 

V 1 / 



The proof is now by induction using the 5^05 AS— construction of L. 
Writing A = Da A^ and B = Da -B* (with Da denoting the obvious subma- 
trix of finite size of the previously infinite matrix Da) we have 



B A 
B A 
\ B A B A J 



\ I A B B \ 
ABA 
A B 
V A/ 
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/ AA AB AB \ 

BA BB + AA AB BB + AA 

BA BB + AA BA + AB 

\ BA BB + AA AB + BA 2BB + 2AA ) 

We have BB + AA = by induction. The result is now correct up to signs 
if AB = -BA. 

We conjugate now this last matrix by the diagonal matrix with diagonal 
entries sq, • • • , S2"+2-i which we can write using blocks of size 2" in the form: 



/ D, 



V 



Z?. 



-^n 



\ 



D^ J 



We get (after simplification using induction and the assumption AB 
-BA) 



/ DaAAD^ 


DaABD^ 





D^ABD^ \ 


D^BAD^ 





-D^ABDa 







-DaBAD^ 







\ D^BAD^ ) 

Supposing the identity AB = —BA we get 

-DaBAD^ = DaABD^ and - D^ABD^ = D^BADa 

which have the correct form (one's on the antidiagonal and zeros everywhere 
else) by induction. 

a 



We have yet to prove that AB 
b 



b a 



-BA. Writing A 
and A = Da A\B = Da B* we have 



b a 



,B 



and 



AB 



BA 



a 
b a 



b 
b a 



b 
b a 



a b 
a 



ab 

ab 

6a 

bd 



which proves AB = —BA since ab = —bd by induction. 



D 
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7 Proof of Theorem II. Ml 



Denote by D^. the diagonal matrix with diagonal coefficients 1, 0, 1,0,1,0,..., 1, 0, 
and by Do = 1— -De the diagonal matrix with coefficients 0, 1, 0, 1, 0, 1, ... , 0, 1,. 
We have thus 1 = Dg + Do and Da = De — Do- The main ingredient for prov- 
ing Theorem 11.31 is the following proposition which is also of independent 
interest. 

Proposition 7.1 We have 

M De L = A + De, M Do L = A + Do 

where A is the strictly lower triangular matrix with coefficients aij = 1 if 
i > j and aij = otherwise. 

Proof. We denote by Oij = (M Df> L)ij respectively /3jj- = (M Do L)ij 
the corresponding coefficients of both products. The proof is by induction 
on i -\- j. Obviously, Oij = Pij = for i < j. The proof for i = j is obvious 
since M and L are lower triangular unipotent matrices. Consider now for 
i > j 





X^ fi + 2k\ (Ak + 1\ 


(recall that 
(mod 2)). 
Since 


-.^ = {'t^) - {'t^) (mod 2) and h, = ci:^}) ^ ci:^}) 


/2i + 2A;\ 
[ Ak )~- 

we have 


/2i + l + 2k\ , / 4A; \ / 4A; + 1 \ / Ak + 1 \ 
~-\ Ak J '''' \2k - 2j) - \2k - 2j) - \2k - (2j - 1)) 




^fli/2\+k\f 2k \ 
""^■"VV 2A: ){k-[{j + l)/2\) 


= 


y. f[i/2\+k\r 2k + 1 \ 

^ \ 2k J\k-\(j + l)/2\y 
fc, fe=L(i+l)/2j (mod 2) ^ ^ ^ LU 1 ;/ J/ 



aii/2\,l{j+i)/2\ if L(i + 1)/2J = (mod 2) , 
P[i/2i,[U+i)/2i if L(J + 1)/2J - 1 (mod 2) . 
We have similarly 



_^^ fi + 2k + l\f Ak + 3 \ 



k 



^ /2[(i-l)/2j +1 + 2A; + 1\ / Ak + Z 
^\ Ak + 2 y V2A; + 1-2[j72J 
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^ n{i-i)/2\+k + i\/' 2k + 1 \ 

fc, fc=[(i-l)/2j (mod 2) ^ / \ U/ J/ 



[{i + 1)/2J +k\f 2k + I 

2k J\k- [j/2\ 



fc, fc=[(j-l)/2j (mod 2) 

"L{i+i)/2j,b72j if L(^ - 1)/2J = (mod 2) , 
Pli^+l)/2^,U/2\ if L(^ - 1)/2J ^ 1 (mod 2) . 
except if [{i — l)/2j = [j/2j = s where 



^ \ 2k + l )[k-s J ~ 

k, k=s (mod 2) ^ ^ ^ ^ 

since only k = s yields a non-zero contribution. This ends the proof. □ 

Proof of assertion (i) in Theorem II. 3I Using the notations of Propo- 
sition [T^ we have 

M DaL = (M De L) - (M Do L) = {A + D^) -{A + Do) = Da 

which shows that Da M Da = L^^ and ends the proof. □ 

Proof of assertion (ii) in Theorem 11.31 We have 



i + j\ (i+ J + 2"+^ 



for i,j < 2"^"^. This shows (together with induction) that the lower right 

[A' \ 

corner of M(2"+^) is given by I „/ ,, ] and has thus the correct form. 

In order to prove the recurrence formula for the lower left corner, we 
remark that B' is by induction the horizontal mirror of A' . We have thus to 
show that ?Tij_|_2n+ij = ?7i2"+i-i-ij or equivalently that 

'2"+i +i + j\ _ /2"+i -l-i + f 



V J \ 2J ■ <°""'^) 



for <i,j < 2"+i. 
Consider 



(2«+i + i + j)(2"+i +i + j-l)--- (2^^+! +i-j + l) 
for < i, J < 2""'"^. Since all terms of the numerator are < 2""'"^ we have 



>"+l +i+j\ _ (2"+i+i+i-2"+^)...(2"+i+i-i+l-2"+2) . , „^ 

2j ) = m^- (mod2j 

(2"+l-i-i)(2"+i-t-i+l)-(2"+i-»+j-l) ^^^^ 2) 



5"+l_,-_i~l('9"+l-i-j-|-n...l'9"+l- 

(2"^^-]-+0 (mod 2) 
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which ends the proof. □ 

5^0S^S-proof of assertion (i) in Theorem [T31 ComDuting L(4) Da M(4) 
(with Da denoting a diagonal matrix with alternating entries 1,-1,1,-1,... 
of the correct size) we have 



/I 



\ 



1 1 
1 1 

Vi 1 1 ly 



/I 



\ /I \ 

1 1 

1 1 1 

-1 / V 1 1 1 / 



Da. 



Writing 



\ 



we have 



B' A' 
A' & A' 
\ & B' A' J 



( ^ 

B A 

{) B A 

\ B A B A ) 



\ 



Da 



I A' 

B' A' 

A' B' A' 

\ B' B' A' ) 



B' A' 
A' B' A' 
\ B' & A' ) 



\ 



/ AA' 

BA' + AB' AA' 

BB' + AA BA' + AB' AA' 

\ 2BA' + 2AB' AA' + BB' BA' + AB' AA' J 

(Da \ 

Da 

BB' + AA' Da 

\ AA + BB' Da J 

by induction. 

Since we have by induction AA' = Da, it will be enough to show that 
BB' = -Da. Writing 



B 



b 
b a 



and B 



we get 



BB' 



b 
b a 



a' b' 
b' 



^ a' 6M 




Kb' J 




bb' 





ba! + ab' 


bb' 



Since ba! + a6' = by the computation above and bb' = —ad' by induction, 
we get the result. □ 
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8 ML and LM 



As before, we denote by L and M the lower triangular unipotent matrices 
with coefficients lij,mij S {0, 1} defined by 



'«j 



T-/) Hod2), m,,^('+.^ 



(mod 2) 



for < i,j. 

Proposition 8.1 The matrix R = M L has coefficients 

ifi<j 



ifi=j 
ifi>j 



Proposition 8.2 The product LM yields a matrix which can recursively he 
constructed by iterating 

\ 



A 
B A 



( ^ 



B A 
2A B A 
\2B 2A B A j 

starting with A = \ and B = 2. 

Assertion (i) of Theorem 11.31 implies of course easily the identities 

(M L)-^ = Da M L Da and (L M)~^ = Da L M Da . 

Proof of Proposition 18.11 This follows immediately from Proposition 
rm since ML = (M D. L) + (M D^ L). □ 

Proof of Proposition 18.21 Using the i3A0i?yl-B— construction of L 
and M we have 

f A \ / A' \ 

B A B' A' 



A' B' A' 
\ B' B' A' ) 



Q B A 

\ B A B A ) 

I AA' 

BA' + AB' AA' 

BB' + AA' BA' + AB' AA' 

\ 2BA' + 2AB' AA' + BB' BA' + AB' AA' J 



\ 



Computing 



AA' 



a 
b a 



and 



BB' 



a 

b' a' 

a' b' 



aa 
ba' + ab' aa' 



b\ f a' b' \ _ f bb' 

b a ) \ b' ) ~ \ba' + ab' bb 

shows AA' = BB' by induction. This finishes the proof by induction. 



D 
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9 Proof of Theorem 11.41 

We prove Theorem 1 1 . 41 bv showing that the n x n submatrix formed by the 
first n rows and columns of the Stieltjes matrix associated with YlT=o ^^ 
has determinant s„ for all n. 

Proof of Theorem 11.41 Let H and L be as in Theorem ll.2l We denote 
by L_ the infinite matrix obtained by deleting the first row of L. Denote 





1 if (i-j) e {-1,0} and c. 



«j 



by C the matrix with coefficients Cij 
otherwise (for < i,j). 

Elementary properties of the Catalan triangle show now that L = L_ C. 
This implies that the finite matrix L-{n) formed by the first n rows and 
columns of L- has determinant 1. Denoting by -Ds_,_i(n) the n x n dia- 
gonal matrix with entries si,. . . ,Sn we see that the finite Stieltjes matrix 
S{n) = Dg{n) {L{n))^^ D s{n) Dg^i{n) L^{n) Ds{n) associated with H has 
determinant s„ G ±1. 

Easy computations show that det{H{2)) = —1 and det{H{2^)) = 1 for 
< k ^ 1. We have also det{H{2^ + a)) = (-1)'^ det{H{2^ - a)) for 
< a < 2^. This shows that det{H{n)) = (— 1)V2J and implies that the 
Stieltjes matrix associated with H is of the form 



/ 1 



Sin) 



1 

d2 1 
-1 ds 1 
-1 d^ 



\ 



We have thus 

Sn = det(5(n)) = d„ det(5(n - 1)) + det(5(n - 2)) = dnSn-i + Sn-2 
for n > 1 which shows 

•Sji Sn—2 



dr, 



Sn-1 



€ {0,±2}, n> 1 



and proves the result. 
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10 The shifted Hankel matrix 

We consider the shifted Hankel matrix H with coefficients hij € {0, 1}, < 
1 if i + J + 2 = 2 for some natural integer k > 1 
and hij = otherwise. We describe the LC/— decomposition of H. The 
associated Stieltjes matrix can be recovered from Theorem 11.11 by setting 
Xi = X2 = xs = . . . = x^ . 



,j given by h, 



*j 
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Define lower triangular matrices L and M with coefficients /j j , rhi j & 
{0, 1} given by k, ^ (^) - (^^^J = (^) (mod 2) and m, = (^W) 
(mod 2). 

It easy to see that one has recursive formulae 

l2i,2j = h,j = ( • _ ■ ) (mod 2), ?2j+i,2j+i = h,j 



and 

'i+j 

2j 



?7i2i,2i = ruij = „ . (mod 2), m2i+i,2i+i = ■m-ij • 



The products L M and M L satisfy analogous recursive formulae. 

In order to state the main result of this section we need also two sign- 
sequences so,si, . . . and to,ti with Si,ti £ {±1}- The sequence Sj is recur- 
sively defined by 

S2i = {-ly and S2i+i = Si 

and is obtained by removing the initial term of the famous paperfolding 
sequence 

1,1,-1,1,1,-1,-1,1,1,1,-1,-1,1,-1,-1,... 

describing the peaks and valleys in a strip of paper which has iteratedly 
been folded (with all foldings executed in a similar way), see e.g. PP, [2] or 

The sequence U is defined by to = 1 and 

^21+1 = ti, tii = { — iyt2i, t4J+2 = t2i ■ 

As always, we denote by Dg, respectively D^ the diagonal matrices with 
diagonal entries sq, si, . . . respectively io,ii, . . .. 

Theorem 10.1 (i) We have 

Di L Ds L* Di = H . 



(ii) We have 



L D-s M = M D's L = D's 



Proof. Let r^ j be the coefficient («, j) of the matrix product Di L Dg U D^. 
Since kj = if i ^ j (mod 2) we have easily that rij = if i ^ j 
(mod 2). For r2i^2j the computation reduces to assertion (i) of Theorem ll.2l 

since l2i,2j = {2i-2j) = Ci-j) ^nd t2j = Sj- 
For r2i+i,2j+i we get 



hi+lS2j +l J2k { 2i+l-2k -l)^^k+l (2j+/-2fc-l) 



i—kJ '^\j~k, 
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which proves assertion (i) by induction on i + j. 

In order to prove assertion (ii), we denote by rij the coefficient {i,j) of 

the product L D^ M. For i ^ j (mod 2) we have kj = ( ^^- ) = rhij = 



i^2i+i ) ~ ^- This shows that rjj- 



if i ^ j (mod 2). For r2i^2j we have 



E 



4i + 2 



2i-2k 



S2k 



2k + 2j + 1 



= E 

k 



2i + l 
i — k 



4j + l 

kik+j\ r 



2j 



if i / i 



by Theorem 11.31 which proves assertion (ii) in this case. 
For r2i+i,2j+i we get 



V V2i + 1 

k 



4i + 4 
2A: 



S2k+1 



2k + 2J + 3 



4j + 3 



^ /2i + 2\ /^A; + j 



2j + l 



= r. 



«j 







S2i+1 



which proves the result by induction on i+ j. 



iii= j 



D 



10.1 SylOSAS-constructions for L and M 

Denote by Lq the infinite strictly lower triangular matrix obtained by adding 
a first row of zeros to L. We define also Mq by adding a first row of zeros 
to M. 

It can be shown that the matrix Lq is obtained by iterating 



A \ 


( ^ 
B 


\ 

A 


B A)^ 





B A 




V^ 


aba) 



where one starts with with 



A A _ / 
B A J ~\l 
Similarly, Mq is obtained by iterating 



A 
B A 



starting with 



A 
B A 



(A \ 

B A 
ABA 
\ B B A ) 





1 
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11 Uniqueness 

Sequences with generating function 

oo 

^efc x^", EiG {±1} 
fc=0 

are bijectively related to the set of all paperfolding sequences and can be 
characerized in terms of Hankel matrices as follows: 

Proposition 11.1 Let so,si,... be a sequence with Si G {±1,0} such that 
det{H(n)), det{H{n)) = ±1 for all n where H{n) and H{n) have coefficients 
hi J = Si^j and hij = Sj+j+i for < i,j < n. Then YlT=0'^k x''^^ = 
X^fcLo ^k x^ for a suitable choice eo, ei, . . . G {±1}- 

Proof. Computing det(i/(n)) (mod 2) determines S2n-2 (mod 2) and 

computing det{H(n)) (mod 2) determines S2n-i (mod 2). This shows 

that ^^0 "^^^ x^^ = Ylk^=o ^ (mod 2) and the signs can of course be 

chosen arbitrarily. D 

I thank J. P. AUouche, M. Mendes France and J. Shallit for helpful re- 
marks and comments. 
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